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Abstract 

Black objects lose their mass and angular momenta by Hawking radiation. It was shown that if 
only scalar fields are emitted, a five-dimensional Myers- Perry black hole evolves toward a state with 
nonvanishing nondimensional rotation parameters, similarly to the case of a four-dimensional Kerr 
black hole. In this paper, we extend this analysis to the case of a black ring. By approximating 
a thin Pomeransky-Sen'kov black ring by a boosted Kerr string, we derive a set of equations that 
determines time evolution of quasistationary evaporation of the black ring. Our result shows that a 
thin black ring evaporates without changing the nondimensional rotational parameter along S . In 
the case of an Emparan-Reall black ring, we solve those equations semianalytically and determine 
its time evolution. 
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I. INTRODUCTION 



In four spacetime dimensions, a stationary, asymptotically flat, vacuum black hole is 
completely characterized by its mass and spin [1]. In particular, the topology of its event 
horizon must be a sphere 2j. By contrast, in five dimensions, in addition to the Myers- 
Perry black hole |3[ which is natural generalization of the four-dimensional Kerr black hole, 
various exact solutions of black objects with nonspherical horizon topologies have been 
found (see J4] for a review). In this paper, we focus attention to the black ring solution with 
the S 1 x S 2 horizon topology. The black ring solution rotating in the direction of S l was 
found by Emparan and Reall j^J. Since a five-dimensional spacetime can have two angular 
momenta, Pomeransky and Sen'kov jsl extended it to black rings with two independent 
rotation parameters (i.e., doubly spinning in directions of S l and S 2 ). 

A black hole is known to evaporate due to the effect of quantum field theory in curved 
spacetime as shown by Hawking [7fl. The rate of mass and angular momentum loss of a 

n n 

Kerr black hole by the Hawking radiation was first studied by Page [8|, |9[ taking account 
of fields with spins 1/2, 1, and 2, and it was shown that a Kerr black hole spins down to a 



nonrotating black hole regardless of its initial state. However, Chambers et al. 



10( (see also 



showed that if only massless scalar fields are taken into account (i.e., in the absence 



of fields with nonzero spin), a four- dimensional Kerr black hole evolves to a state with the 
nonvanishing nondimensional rotation parameter, a/M ~ 0.555. This analysis was extended 
to five-dimensional Myers-Perry black holes by Nomura et al 12] . They showed that any 
such black hole with nonzero rotation parameters a and b evolves toward an asymptotic 
state with a/M 1 / 2 = 6/M 1 / 2 ~ 0.1975(8/3tt) 1 / 2 . Here, this value is independent of the 
initial values of a and b. 

ack ring. Although the Hawking 
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16) , so far, the time evolution 



It is interesting to extend these studies to the case of a b 
radiation of black rings has been studied in various context 
of a black ring has not been studied. The difficulty in this study is that the method of mode 
decomposition of the Klein-Gordon field in this spacetime is not known since separation 
of variables has not been done, and therefore, two dimensional numerical calculations of 
eigenfunctions are required. In order to avoid this difficulty, we take the thin limit and 
approximate a black ring as a boosted black string. In this approximation, the variables can 
be separated and we have well defined modes. 
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Using this approximation, we give a formulation to study the evolution of a thin 
Pomeransky-Sen'kov black ring by the Hawking radiation, and discuss general features that 
do not depend on details of greybody factors. Then, we derive a semianalytic solution in the 
case of an Emparan-Reall black ring and determine its time evolution. Here, the solution is 
semianalytic in the sense that one parameter that is related the greybody factors has to be 
evaluated numerically, although the evolution is expressed by analytic formulas. 

This paper is organized as follows. In Sec. II, the black ring solution is reviewed and 
its boosted Kerr string limit is shown. In Sec. Ill, we derive the equations that determine 
the emission rates of mass and angular momenta of a black ring via Hawking radiation. In 
Sec. IV, the evolution of a black ring is discussed. Sec. V is devoted to a summary. In 
Appendix A, we check the validity of our numerical result by studying the DeWitt approx- 
imation. To simplify the notation, the natural units h = c = G = kB = l are used. 



II. BLACK RING 



In this section, we review basic properties of a black ring and show its boosted Kerr string 
limit. This limit was discussed in the more general cases of an unbalanced Pomeransky- 

n 

Sen'kov black ring in Ref. |17| . 

The metric of the Pomeransky-Sen'kov black ring is {f| 



ds* 



+ 



H(x,y) y ' "' H(y,x) r ~H(y,x) 
2R 2 H(x,y) ( dx 2 dy 2 \ 



H(y,x) 



{x _ y f (i _ v f \ G {x) G(y)J 
where the 1-form Q is 



Q = - 



2RX\ (1 + vY-X 2 



H(y,x 
1 + y 



(l — x 2 ) yy^udcj) 



1-X + v 



{1 + A - v + x 2 yu(l - X - v) + 2vx (1 - y)} dtj) 



and the functions G, H, J and F are 



(2) 



G(x) = (1-x 2 ) (l + Xx + vx 2 ) 



(3) 



H(x, y) = 1 + A 2 - v 2 + 2Az/ (l - x 2 ) y + 2xX (l - y 2 v 2 ) + x 2 y 2 v (l - A 2 - v 2 ) , (4) 
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2R 2 (1 — X 2 ) (1 — V 2 ) Xjv r 9 9 , , , / o 9\ T / v 

J(x, y) = \ !± y -t^ + \ —v +2(x + y)\v — xyv (l — A — v 2 )\ , (5) 

(x - y) (1-p) 

b [x, y) = g 

{x-y) (1 - v) 

G(x) (1 - y 2 ) [{(1 - vf - A 2 } (1 + v) + y\ (l - A 2 + 2v - 3u 2 )] 

(6) 

+G{y) [2\ 2 + x\ {(1 - v) 2 + A 2 } + x 2 {(1 - v) 2 - A 2 } (1 + v) 

+ x 3 A (1 - A 2 - 3u 2 + 2v 3 ) - x\l - v)v (-1 + A 2 + v 2 ) ] 

Here, we follow the notation of Ref. [6| except that we choose the signature (— , +, +, +, +) 
for the metric, exchange <fi and ip, and use R instead of k. The coordinate ranges are 
— oo < t < +oo, < (ft, ip < 2n, — 1 < x < 1 and — oo < y < —1. R is a parameter of 
dimension of length, which determines the characteristic scale. A and u are dimensionless 
parameters satisfying < v < 1 and 2\fi> < A < 1 + which determine two nondimensional 
rotational parameters. The regular event horizon exists at y = yh, where 



_A+VA^4^ 

Vh = 2v • (7) 

The solution is asymptotically flat and spacelike infinity is located at x — y = — 1. The 

coordinates (x,(f>) parametrize the two-sphere S 2 and ip parametrizes the circle S 1 . One 
recovers the Emparan-Reall black ring by setting v = 0, and the line v = A 2 /4 represents 
the sequence of an extremal black ring (see Fig. [1]) . 
The mass and angular momenta are 

_ 3nR 2 X _ 4n&\^y/(l + v) 2 - A 2 

M ~ 1-A + z/ J +~ {l-v) 2 {l-\ + v) ' (8) 



_ 27ri? 3 A(l + A - 6u + \v + v 2 ) ygj + v) 2 - A 2 
*~ (l-z/) 2 (l-A + z/) 2 ' ' > 

The angular velocities, the area, and the surface gravity of the horizon are written as 



A( i + i/ )_(i_, )v ^; / 1 + 1/-A i / 1 + I/-A 

n * = ^ VTT^Ta' ^ = 2^VTT^Ta' (10) 



_ 327r 2 i? 3 A(l + z/ + A) fa; 1 - - z^y/A 2 ^!^ 

(l-^fa- 1 -?/,) ' ^ 4i?A(l + z, + A) • 

We consider a thin ring limit A — > where the ratio of the S 2 radius to the S 1 radius 
becomes very small. Here, we have to take care of the fact that this limit depends on the 
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FIG. 1: The parameter space (A, v) of the Pomeransky-Sen'kov solution. A and v can take values 
in the region surrounded by solid lines. The line v = A — 1 is the Myers- Perry black hole limit, 
the line v = is the Emparan-Reall black ring, the line v = A 2 /4 is the extremal limit. The point 
A = v = corresponds to the boosted Kerr string limit. The dashed line v = a 2 A 2 /4 is the path 
to the point A = v = with a fixed a*. 



path to the point A = v = 0. For example, taking limit A — > on the line v = gives the 
boosted Schwarzschild string, while taking limit A — > on the extremal line ^ = A 2 /4 should 
result in a state different from the boosted Schwarzschild string. Therefore, A = v = is a 
degenerate point, and in order to resolve this degeneracy, we introduce a new parameter a* 
as 

v = -al\\ (12) 



and consider the limit A — > on the line of a fixed a* (see Fig. [T]). Also, in order to obtain 
a well-defined limit, we introduce 



M K = —XR 



V2 



(13) 



and fix Mr in taking this limit. 

We introduce the new coordinates r, z and 9 by 

V2R 



V2R'' 



x = cos 6, 



(14) 



and collect the leading-order term of each metric component with respect to A. Then, the 
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black ring solution is reduced to the so-called boosted Kerr string solution 
^ = _ , 1 _ 2M K rco^a\ ^ + 2M K r^2a ^ + / + 2MggnhV | ^ 



P 2 , 2 2 ,/,2 ( r2 + fl2 ) 2 - Aa 2 sin 2 9 . ., .. , , 

+ ^—dr 2 + p 2 d6 2 + - - — sin 2 u 1 - 

A p 2 



AMkt cosher . 2 4Mxrsmhcr . 2 

asm adtdep asm odzdep, (15) 

p 2 p 2 

where p 2 = r 2 + a 2 cos 2 9 and A = r 2 — 2Mkt + a 2 , and a is defined by a := MkCl*- Since 
is the mass density and a is the Kerr parameter, a* represents the nondimensional rotation 
parameter along S 2 direction, a := arctanh(l/\/2) is the boost parameter. Although the 
boost parameter can take any value for general boosted Kerr strings, it is restricted to this 
value for the thin limit of a black ring. The event horizon is located at r = r + where 



r± := M K ± v/M£ - a 2 . 

In studying the evaporation of a black ring, we use this boosted Kerr string solution in 
the following sense. We consider the situation where A is very small, and do not take the 
exact limit. Then, in the neighborhood of the black ring, the spacetime metric can be well 
approximated by the boosted Kerr string solution. For this reason, the value of R is not 
infinite in our analysis although it is very large compared to Mr. 

In this thin-limit approximation, the physical quantities in Eqs. (18])- (lll|) are expressed 
in terms of R, Mk and a (or a*) as: 

M ~ 3V2irRM K , ~ 2V2ttR 2 M k , ~ Aira^RM^. (16) 

A H ~16n 2 r + M K R, k — ~tvt~~~ t • (18) 

4Mxr + cosh a 

The inverse relations of Eq. f TTB]) are 

K ~ 2M' K ~ 9tt J 4 , ' * " 4 M3 ' liyj 

and A is expressed as 

V2M K 4 M 3 

A ~ ~ (20) 

R 27tt J 2 1 ; 

Because A -1 / 2 is proportional to the angular momentum J$ normalized by mass M, A can 
be interpreted as the nondimensional rotation parameter along S l . 



Note that in Eq. ffTTl) agrees with the angular velocity defined by the Killing generator 
£ of the horizon of the boosted Kerr string 1 , 



£ = dt + + vd z 



(21) 



with 



V = tanh er, 



(22) 



and k in Eq. (|18p is in agreement with the surface gravity of the horizon calculated with £. 
For a later convenience, it is useful to compare and k with the angular velocity and the 
surface gravity of the horizon of an unboosted Kerr string. In the following, the quantities 
in the unboosted system is indicated by prime ( ' ). In the unboosted system, the Killing 
generator of the horizon is £' = d t > + ^<9</>, and Q!i and k' calculated from £' are 



There are differences in the quantities of the boosted and unboosted systems by a factor of 
1/ cosher. This is understood as the effect of time delay in the Lorentz boost. 

III. FORMULATION 

In this section, we formulate the emission rates of mass and angular momenta of a black 
ring via Hawking radiation. 

A. Emission rate 

In this subsection, we derive general formulas to determine the evolution of a thin 
Pomeransky-Sen'kov black ring by the Hawking radiation utilizing the black string approx- 
imation. The evolution of a scalar field is governed by the Klein-Gordon equation in curved 
spacetime 




a 



r + — r_ 
AM K r + 



(23) 



2M K r + ' 



K = 




(24) 



where g is the determinant of its metric. 



Our expression of fi^ does not agree with that of Ref. [19| because the definition is different 
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To quantize the field, we need to expand it in terms of eigenmodes for $, which can be 
written in the black ring background as 

<f> = e- i " t e im(l 'e in ^(x,y), (25) 

where u, m and n are the eigenvalues for the Killing vector fields d t , and d^, respectively. 
By inserting this expression into Eq. (1241) . we obtain a second-order elliptic equation for 
\I/(x, y) in the (x,y) plane. This equation has a discrete series of regular solutions labeled 
by an integer I. In the Schwarzschild string limit with = 0, this series of solutions 
become proportional to the associate Legendre functions P™(x). Thus, the mode functions 
are labeled by the four parameters (u, £, m, n) in which £, m and n take integer values. 

In the case of a nonrotating black hole, the expected number of particles emitted per unit 
time is proportional to black body radiation: 

where u s is the energy of a scalar field in the background of the nonrotating black hole and 
T s is the temperature of the horizon. Here, the temperature is defined as T s = k s /2it, where 
k s is the surface gravity of its horizon. The temperatures appearing later are defined in a 
similar way. 

In the rotating case, we have to replace u s by the energy of the mode with respect to 
the null geodesic generator of the black hole horizons because the mode function behaves 
as exp(— iu*u±) in the coordinates that are regular around the black hole horizon, where 
u± is the advanced time/retarded time around the horizon. In general, this null geodesic 
generator £ can be written as £ = Qt+J^ . Qjd^ in terms of the time translation killing vector 
d t and the rotational Killing vectors . Hence, for the mode oc exp(— uot + i ^ . m^) 
is expressed as 

= to — wPVtj. (27) 

3 

Hence, the expected number of particles emitted per unit time from the black ring is given 
by 

(-^br) = e(w _ n n^ )/T B R ( 28 ) 

where Tbr is the temperature of the horizon and r^^(u;) is the greybody factor, which is 
identical to the absorption probability of the incoming wave of the corresponding modes. 
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This leads to the emission rates of the total mass and angular momenta, which give changes 
of the black ring mass M and angular momenta and as 



d 
di 



■I 



-E 

£,m,n 



dco(N) 



n 



(29) 



where E^ means summing up over all modes. Note that in this expression, it is difficult 
to estimate the greybody factor generally because no one has succeeded in separating the 
coordinates x and y, and two-dimensional numerical calculation of eigenfunctions is required. 

Here, we consider the situation where the mode functions can be approximately evaluated: 
a black string limit discussed above. For the boosted black string (fl5|) . we can separate the 
wave equation, and therefore, we approximate the evolution of a scalar field in a black ring 
spacetime by that in a boosted Kerr string spacetime. In this situation, the scalar field can 
be separated as 



-iuit 



R(r)e 



-ikz imd> nm 



^st{o), 



(30) 



where S™(9) is the spheroidal harmonic function. From the coordinate transformation ffT4j) . 
the relation between n of a black ring and k of a boosted black string is n = 2fci? tanner. 
The expected number of particles emitted per unit time is given by 



(AT) 



(31) 



e (uj-kV-m,n^)/T _ J_> 

where T = k/2tt is the temperature of the horizon with k in Eq. ( TT8|) . and V is the linear 
velocity of the horizon along the string direction, Eq. f[2"2"j) . Ti mn (u) is the greybody factor 
of the boosted Kerr string spacetime. The emission rates (l2"9l is evaluated by using (N) 
instead of (A^br). 



B. Simplification 

In the following, we simplify the expression of the emission rates (1291) . We normalize all 
quantities by the mass density Mk, 

Co = M K u, k = M K k, % = MkQ^, f = M K T. (32) 

The normalized quantities Q'^ and T" in the unboosted system are introduced in the same 
manner. It is important to point out that the normalized angular velocity and the 
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normalized temperature T depend only on a*. We change the order of summations over i 
and m by 

oo I oo oo 

££=££. (33) 

^=0 m=—t m=—oo £=\m\ 

and define 

oo 

S (m) M) := ^ r lmn . (34) 

£=\m\ 

Then, the emission rates can be written as 



d 
dt 



J.,j, 



2-kM k 



£ £ 



m=— oo n=— oo 



doo- 



' Cj/M k ^ 



fc| 



^Q-kV-mil^/T 



\ 



n 
m 



(35) 



/ 



Here, the lower limit of the integral is \k\ because the mode with its energy oo < \k\ is 
gravitationally bounded and does not escape to infinity. Because the number of the mode 
n per unit k is very large, 0(1/ X), the summation over n can be replaced by the integral: 

2R tanh a 



£ 



— > 



dn 



M, 



K 



dk. 



(36) 



Because the integrand of Eq. (135 j) is zero on the boundary of the integration domain, the 
difference of the integral from the original discrete summation is 0(X 2 ) and thus negligible. 
Rewriting Mk and R by M, Jf and using Eq. f JT9|) . we obtain 

Cog {m \Co', k') 



dM 


3 ? 7T 2 4 


dt 


2 3 M 7 


dJ^p 


3 8 vr 2 4 


dt 


2 3 ^2M 8 


dJ^ 


3 5 7T 4 


dt 


2 2 V2~M 5 



£ 

m=— oo ' 
oo 

£ 

ra=— o 
oo 

£ 

m=— oo ' 



m=— oo 



(ifc 



f/A' 



(id) 



(id) 



kg( m \co',k') 
mg^(Co',k') 



(37) 

(38) 
(39) 



The integral in each formula can be further simplified if we perform transformation from 
(oo, k) to (oo', k ), where 



oo = Co' cosh er + k' sinh a, k = oo' sinh a + k cosh a. 
Substituting these formula with cosher = y2 and sinh a = 1, and defining 



/{ m) (a*) 



dk' 



doo 



r , Co'g^(Co',k') 



\k>\ 



3 (u'-mfly/T' 



(40) 



(41) 
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we obtain 



/ Ar,i u'-ran'J/T' -i ' v ' 



M (it v y M 8 ' 



GM-rb*> ( 45 ) 



with 



J 6 dt v 'J^M 5 



F ^) : =|^f E J i M K), (46) 



m=— oo 



G ^ ■= £k E ^V*)- (47) 

v m=— oo 

Note that we used the fact that the integral of terms proportional to k vanish because such 
terms are odd functions of k . From Eqs. fj43|) - fj45l) and Eq. ffT9]) . the equation of a* is 
derived as 

1 ''" 3H(«.)4. (48) 



where 

9ti G(a*) „, . . 

if(a«) := — -AJ^_f(o»). (49) 

Therefore, the time evolution of a black ring in the Hawking radiation is determined by 
the equations for M, and a*, that is, Eqs. ( H3|) . (jUJ) and (|48l) . The remaining work 

is to calculate the greybody factors and obtain F(a*) and H(a*) of Eqs. fj46]) and f l49|) 
numerically. 



C. Greybody factor 



In the following, we discuss the greybody factors for massless scalar fields in a boosted 
Kerr string spacetime. Substituting the ansatz ( 1301 into the Klein-Gordon equation (I24|) 
in the background of the boosted Kerr string (IToT) . we get the following angular and radial 
wave equations for S™(9) and R(r) [l9j , 







sin 9 



d e {sm9d e S™) + 



(u 2 - k 2 ) cos 2 9 —j- + Xtm 

sin' 1 9 



SI 



(50) 
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= Ad r (Ad r R) - A [k 2 r 2 + a 2 u 2 - 2wma cosher + \e m ] R 

+ \u (r 2 + a 2 ) — ma cosh <r] 2 + 2M^r (r 2 + a 2 ) cosh 2 a (u — k tan a) 2 

- 2M K r (r 2 + a 2 ) w 2 - m 2 a 2 sinh 2 a + 4kmaM K r sinh a R, (51) 

where A£ m is the separation constant which is determined as an eigenvalue of (|50|) . For 
small a 2 (oo 2 — k 2 ), the eigenvalues associated with the spheroidal wave functions S™ are 
\ lm = t (£ + 1) + O (a 2 (oo 2 - k 2 )) Q. 

As we mentioned in Sec. Ill A, the greybody factor is calculated as the absorption 
probability of the incoming waves of the corresponding mode. With Eq. (jlCTl) . the tortoise 

coordinate r* and a new wave function u defined as 

r 2 _|_ a 2 i 
dr* = — x dr, R = / o = u, (52) 



the radial wave equation f[5"Tj) can be rewritten as the equation of Schrodinger type: 

i 2 

— + u' 2 -V(r) u = 0, (53) 

■ r * 

where V(r) is the effective potential 

A {2M K r(r 2 - 2a 2 ) + a 2 (r 2 + a 2 )} A(u/ 2 a 2 + A £m + k' 2 r 2 ) + Amu'M K ar - m 2 a 2 



V(r) 



(r 2 + a 2 ) 4 (r 2 + a 2 ) 2 

(54) 



As the boundary condition, we impose the ingoing condition at the horizon. Then, u behaves 
as 

{e~* w + r * at r — y r+, 

(55) 
A m e~ iuJ ™ r * + A out e iu] '~> r * at r ^ oo. 

Here cu^ := u' — mQ'^ and u'^ := a/ cj /2 — k' 2 . The greybody factor is written as 

r = 1 - \A out /A m \ 2 } (56) 

which has to be evaluated numerically. 

In this paper, we perform numerical calculations of the greybody factors in the case of 
the Emparan-Reall black ring, i.e., a = 0. In this case, S™ becomes the spherical harmonic 
function and the calculation becomes much simpler compared to the case a ^ 0. We devel- 
oped a code to calculate the greybody factors in this situation. Here we show the example 
of the behavior of the greybody factors for the modes of three lowest I numbers in Fig. |2J 
Our result is in g ood agreement with the analytic approximate formula for low-frequency 



waves in Ref. 
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FIG. 2: Greybody factors for the modes k = m = 0, £ = 0, 1, 2 (from left to right). 
IV. EVOLUTION BY EVAPORATION 

In this section, we discuss general features of time evolution of the evaporating 
Pomeransky-Sen'kov black ring. Then, we focus attention to the case of the Emparan- 
Reall black ring, and derive semianalitic solutions of the time evolution using the numerical 
results of the greybody factors. 

A. Evolution of Pomeransky-Sen'kov black rings 

Here, we discuss the general feature that does not depend on the details of the greybody 
factor. From Eqs. ( H3|) and fPB]) . the following relation can be immediately found: 



Therefore, a Pomeransky-Sen'kov black ring evaporates without changing the initial value 
of the nondimensional rotation parameter along S 1 . 

Next, let us assume a* > and consider how to solve a*(t) and M(t). Eliminating M 
and J$ from Eqs. (jl3"j) and (HBl) using Eq. fl57|l . we obtain 




(57) 



From Eq. (I20p . this means that the parameter A is unchanged, 



A(t) = A(0). 



(58) 




(59) 
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In principle, this equation can be solved at least numerically once if (a*) and F(a*) are 
generated. Then, from Eqs. (1431) and (T4*8i) . the time evolution of M(t) is formally given as 



M(t) = M(0) exp 



F(a, 



o»(0) 



-da* 



(60) 



We point out that the behavior of H(a*) is crucial for the evolution of a*. This function 
is analogous to h(a*) of Ref. lOj where evolution of a four-dimensional Kerr black hole was 
investigated: The value of a* increases (decreases) if if (a*) is negative (positive). If H(a*) 



crosses zero from negative to positive at some value a* = a* similarly to Fig. 3 of Ref. 10]. 

(c) 

the black ring inevitably evolves to the state with a* . Therefore, numerical calculation of 
H(a*) is very interesting and we are planning to study this issue. In the next subsection, as 
the first step, we study what happens in the case of an Emparan-Reall black ring. 



J^t) = J40)[l-AF(0)^-t) . (62) 



B. Evolution of Emparan-Reall black rings 

In the case of the Emparan-Reall black ring, </</,(£) = and a* = and Eqs. fT43]) and fj4*8|) 
can be solved analytically: 

M(t) = M(0) (l - 4F(0)^|^f) V2 , (61) 

^(Q) 4 x3/4 
M(oy 

Therefore, our remaining task is to determine F(0) numerically. 

As discussed in Sec. IIIIBt basically we have to calculate Eqs. ( 14T|) . and fT46|) . and in 
those equations, the summation over m was taken at the last step. But in the case of a* = 
considered here, it is better to take summation with respect to m in advance because the 
integrand does not depend on m. For this reason, we calculate summation of the greybody 
factors as 

oo l oo 

g(u', k') = J2 J2 T 'l™ = + ( 63 ) 

£=0 m=-l 1=0 

Then, -F(O) is given by 

m ^ h , /i: .r^r^% (64) 



2 3 y/2 ' J-oo Jw\ e*'/ fl - 1 

The integration of I\ was proceeded with the Simpson's method. The domain of integration 
was made finite by discarding the region where the integrand is smaller than 10 -16 . We took 
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summation with respect to £ up to £ = 10, because the greybody factors for £ > 10 turned 
out to be negligible. This is because the potential wall for £ > 10 is so high that it reflects 
almost of all waves. In this manner, F(0) is determined as 

F(0) ~ 0.239. (65) 



As a check, we compute F(0) using the DeWitt approximation 22j in Appendix |A] The two 
results agree well, and therefore, our numerical result is reliable. 

As shown in Eq. fl58l) . the nondimensional parameter A is held fixed throughout the 
evolution. Because a* = 0, the constancy of the nondimensional parameter also indicates 
that the Emparan-Reall black ring evaporates keeping similarity to its initial shape: The 
black ring at any time can be obtained by uniformly scaling its initial configuration. The 
scaling factor can be found by deriving the time evolution of the ring radius R as 

ilrH^^r (66) 

The lifetime T of a thin black ring with mass M is 

/27ttA\ 2 f M x 2 



T *{— ){w p ) t - < 67 > 

where M p and t p are the Planck mass and Planck time, respectively. The time scale is 
proportional to M 2 , and this dependence on M is same as that of the five- dimensional 
Schwarzschild black hole. However, because of the prefactor (277rA/4) 2 , the life time of a 
thin black ring with A < 10~ 2 is much shorter than the five- dimensional Schwarzschild black 
hole with the same mass. 



V. SUMMARY 

In this paper, we have studied the time evolution of evaporation of a thin black ring 
under the assumption that only massless scalar fields are emitted in the Hawking radiation. 
In order to separate the Klein-Gordon equation in the background of a black ring metric, we 
have considered the thin-limit approximation, where the black ring metric is approximated 
by the boosted Kerr string metric. Then, we have given a set of equations, Eqs. (H3"|) . flS]) 
and (JUJ), that determines the quasistationary evaporation of a Pomeransky-Sen'kov black 
ring. In this setup, a black ring evaporates without changing the nondimensional rotation 
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parameter A along the S 1 direction. Also, we have analytically solved these equations in the 
case of an Emparan-Reall black ring, Eqs. (IBTj) and (162|) . with the factor F(0) ~ 0.239 that 
has been determined by numerical calculation of the greybody factor. In the evaporation, 
the shape of the Emparan-Reall black ring keeps similarity to its initial configuration. 

As a future work, we plan to study the time evolution of a Pomeransky-Sen'kov black ring 
with nonvanishing rotational parameter a* along the S 2 direction. For this purpose, we have 
to calculate the functions F(a*) and if (a*) of Eqs. (146]) and (j4~9]) . and therefore, developing 
the code for computing the greybody factors of the boosted Kerr string is required. 
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Appendix A: DeWitt approximation 

In order to check the validity of our numerical result (1651) of F(0) in the case of the 
Emparan-Reall black ring, we compute this value in an approximate way. As the method 
of approximation, we adopt the DeWitt approximation [22] that was originally developed 
to evaluate the contribution of the greybody factor to the evaporation of a Schwarzschild 
black hole (see p. 394 of Ref. [23y for a review). In that study, the greybody factor was 
obtained by appropriately reinterpreting the capture condition of null geodesies. Although 
this approximation holds only for high-frequency regime in a strict sense, it gives a rather 
good result. In fact, the difference of the DeWitt approximation from the numerical result is 
~ 6%. Compare the formula of the mass-loss rate by the DeWitt approximation (Eq. (146) 



of Ref. 22j ) and the numerical value reported in Ref. [lOj . 

In the spacetime of a five-dimensional Schwarzschild string, a massless particle with 
momentum in the z direction effectively behaves as a massive particle in a four- dimensional 
Schwarzschild spacetime. Therefore, as the first step, we study timelike geodesies in a four- 
dimensional Schwarzschild spacetime and derive the capture condition. Then, we translate 
it to the greybody factor for a massless scalar field in a spacetime of a Schwarzschild string. 
Using this result, we derive the value of F(0) in the DeWitt approximation by performing 
the integration and summation in Eqs. (134)) . (14TI) . and (146]) . 
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The metric of a four-dimensional Schwarzschild spacetime is given by 

dr^ 

ds 2 = -f(r)dt 2 + -77-t + r 2 (d9 2 + sin 2 6d<f) 2 ), (Al) 
fi r ) 

m = i - (A2) 

r 

The geodesic motion of a massive particle in the equatorial plane is governed by the following 
equations: 

f(r)t = e, (A3) 

r 2 = I, (A4) 

-f(r)i 2 + ^- + r 2 <p 2 = -l. (A5) 

Here, e and Z indicate the energy and angular momentum per unit mass of the test particle, 
and dot ( ' ) denotes the derivative with respect to particle's proper time r. Substituting 
Eqs. ((3SI) and (TA4]) into Eq. f|A~5]) . we obtain 

f 2 + V(r) = e 2 , (A6) 

where 

V W = ^ + 1 Vl-^V (AT) 



Let us consider the situation where a test particle exists in the neighborhood of the horizon 
and moves toward outside (i.e. f > 0). Denoting the peak value of V(r) as Vp ea k, the 
particle escapes to infinity if the condition Vp ea k < e 2 is satisfied. Conversely, a particle 
with Vp ea k > e 2 is reflected back to the black hole by the centrifugal barrier. After some 
calculation, the condition Vp ea k = e 2 is shown to be equivalent to 

I = e^/F\ejM K , (A8) 

where 

_ 27e 4 - 36e 2 + 8 + e(9e 2 - 8p 
{> 2e 2 (e 2 -l) " 1 ' 



Therefore, a particle escapes to infinity if I < e^J F(e)M K) and it is reflected back to the 
black hole if I > e^J F(e)Mx- This condition is equivalent to the one obtained in Ref. 24]. 
Here, yjF{e) varies from 4 to 3v^3 as e is increased from unity to infinity. 

We use this result in order to approximate the greybody factor in the particle emis- 
sion from the Schwarzschild string. Here, we choose the unboosted frame, and as done in 
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Sec. IIIIB| the quantities in this frame are indicated by prime ( ' ). Consider the emission 
of a quantum particle with mass k', angular frequency u', and angular quantum number £. 
Replacing as I — >• i/k' and e — > uo'/k' in the capture condition derived above, the particle 
with i < u'^F{u'/k')M K escapes to infinity and that with i > u'yj F(uo' / k')Mx falls back 
to the horizon. Therefore, the greybody factor is approximated as 



Cb'dF{Cj'/k')-l), (A10) 



where 6{u) denotes the Heaviside step function, and we introduced u' = M K uo' and k! = 
M^k' in the same manner as Eq. (13"2|) . Note that in the massless particle limit, k' /Co' — > 0, 
Eq. (lAlOp becomes r^ m „ w 6(3y/3u' - £), and this agrees with the formula in the original 



DeWitt approximation 22]. 



Now, we evaluate the value of F(0). The computation can be done with the formula 
given in Sec. IIVBI The quantity g(Cu', k') in Eq. (|63|) is 

oo 

g(Q', k') = ^(21 + l)6(Co'^F(uj>/k') - £) « Co' 2 F{u'/k r ). (All) 

1=0 

Here, the summation over i was changed by integration as done by DeWitt. Then, F(0) can 
be calculated by substituting this formula into Eq. (jMj) . It is convenient to introduce new 
variables (x,y) by x = k' /Co' and y = Co', and in these variables, analytic integration can be 
proceeded as 



C(5) 



88 + 33 V2 arcsin f ) - 3 log 3 



V 3 



(A12) 



4096vr 5 

Substituting this result into Eq. ( IMj) . we obtain 

F(0) w 0.224. (A13) 

This value is fairly close to our numerical value in Eq. (j65|) : Similarly to the original De- 
Witt approximation for the Schwarzschild black hole, the approximate value is about 6% 
smaller compared to the numerical value. Therefore, the result of the DeWitt approximation 
supports the correctness of our numerical calculation. 
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